Interlayer Aharonov-Bohm interference in tilted magnetic fields in 
quasi-one-dimensional layered materials 
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Different types of angular magnetoresistance oscillations in quasi-one-dimensional layered materi- 
als, such as organic conductors (TMTSF^X, are explained in terms of Aharonov-Bohm interference 
in interlayer electron tunneling. A two-parameter pattern of oscillations for generic orientations of 
a magnetic held is visualized and compared with the experimental data. Connections with angular 
magnetoresistance oscillations in other layered materials are discussed. 

PACS numbers: 74.70.Kn, 72.15.Gd, 73. 21. Ac. 



5h 
43 



c3 

a 

o 

(N 
> 

m 
o 

On 
O 

in 
o 
-i— > 



T3 

c 

o 
o 



X 



Angular magnetoresistance oscillations (AMRO), 
where resistivity oscillates as a function of the magnetic 
field orientation, were originally discovered P, m the 
quasi-two-dimensional (Q2D) organic conductors of the 
(BEDT-TTF) 2 X family 0. AMRO are distinct from 
the Shubnikov-de Haas and de Haas-van Alphen oscil- 
lations and are now widely used for direct mapping of 
the Fermi surfaces of layered materials 0, EJ- AMRO 
have been observed not only in many organic conduc- 
tors, but also in intercalated graphite 0, Sr 2 Ru04_ 
Tl 2 Ba 2 Cu0 6 8, 9], and GaAs/AlGaAs superlattices 

Early theories of AMRO H [H E3 

were formulated in 

terms of semiclassical electron trajectories on a cylindri- 
cal 3D Fermi surface. Then it was realized that AMRO 
can exist already for two layers 0,0,0], and they rep- 
resent an Aharonov-Bohm interference effect in interlayer 
tunneling [J^. Some experimental evidence for AMRO 
in semiconducting bilayers has been found 0] , but more 
systematic measurements are necessary. 

AMRO were also found in quasi-one-dimensional 
(Q1D) organic conductors with open Fermi surfaces, such 
as (TMTSF) 2 X 3]. These materials consist of parallel 
chains along the x axis, which form layers with the in- 
terlayer spacing d along the z axis and the interchain 
spacing b along the y axis, as shown in Fig. Origi- 
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FIG. 1: (a) Geometry of electron tunneling between two Q1D 
layers. The sinusoidal line represents the in-plane electron 
trajectory; tt and t c are the amplitudes of interchain tun- 
neling, (b) The Fermi surfaces of the two layers, shifted by 
the vector q in Eq. The shaded areas Si and 52 lead to 
interference oscillation in the presence of B z . 



nally, three different AMRO were discovered in the Q1D 
conductors: the Lebed magic angles 0, 0, H3, El for 
a magnetic field rotation in the (y, z) plane, the Danner, 
Rang, and Chaikin (DKC) oscillations in the (x, z) plane 
23, pj, and the third angular effect in the (x, y) plane 
3SE1- Then Lee and Naughton H3 found com- 
binations of all three effects for generic magnetic field 
rotations. It became clear that all types of AMRO in 
Q1D conductors have a common origin and should be 
explained by a single unified theory. 

In (TMTSF) 2 X, the in-plane tunneling amplitude be- 
tween the chains, tj ~ 250 K [22I I23I ]. is much greater 
than the inter-plane tunneling amplitude t c ~ 10 K [jj]. 
Thus, we can treat interlayer tunneling as a perturbation 
and study the interlayer conductivity a c between just two 
layers for a tilted magnetic field B = (B ^,, B y , B r ), as 
shown in Fig. QJ,. This bilayer approach [lj, |l5j only 
assumes phase memory of interlayer tunneling within a 
decoherence time r and does not require a well-defined 
momentum k z and a coherent 3D Fermi surface. It gives 
a simple and transparent inter pret ation of the most gen- 
eral Lee- Naughton oscillations |2j| in terms of Aharonov- 
Bohm interference in interlayer tunneling. The results 
are equivalent to other approaches based on the classical 
Boltzmann equati on Eill26ll27ll2cil29] | and the quantum 
Kubo formula [3fl l3ll l32j|. We calculate a contour plots 
of a c as a function of two ratios B x j B z and B y /B z for 
models with one or several interlayer tunneling ampli- 
tudes This type of visualization clearly reveals 
agreement and disagreement between theory and exper- 
iment and allows to determine the electronic parameters 
of the Q1D materials. The results can be also applied 
to Q1D semiconducting bilayers consisting of quantum 
wires induced by an array parallel gates, as shown in 
Fig.HJi. 

Let us consider tunneling between two Q1D layers. 
The in-plane electron dispersion is 



s(k x ,k y ) = ±v F k x - 2t b cos(kyb/h) 1 



(1) 



where energy e is measured from the Fermi energy, ±vp 
are the Fermi velocities on the opposite sheets of the open 
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Fermi surface, k — [k x ,k y ) is the in-plane momentum, 
and k x is measured from the Fermi momentum. The 
interlayer tunneling is described by the Hamiltonian 

H±=t c [ ^(r)^ 1 (r)e^ (r) d 2 r +H.c, (2) 



<t>(r) = —A z (r), A z (r) = B x y-B v x, (3) 

where r = (x,y), c is the speed of light, e is the elec- 
tron charge, A z is the vector potential, and ipi 2 are the 
electron destruction operators in the layers 1 and 2. The 
gauge phase <j>{r) is due to the in-plane magnetic field. 

We treat the in-plane electron motion quasiclassically. 
For B z ^ , el ectrons move in time t along sinusoidal 
trajectories |34|, as shown in Fig.QJi, 



x(t) = x ± v F t, y(t) =y ± 



f 2t b c 
\ev F B z 



cos(w c i). (4) 



Instead of the magnetic field components (B x , B v , B z ), it 
is convenient to introduce the variables oj c , B' x and B' y 
defined by the following relations: 



ebvpB z 
he 



BL = 



B x 2i h d 
B z Tivf 



, B», d 

r>f y 

By -^b- 



(5) 



The cyclotron frequency lo c is simply proportional to B z , 
whereas the dimensionless variables B' x and B' y are pro- 
portional to the ratios of the magnetic field components 
B x /B z = cos<ptan# and B y /B z = sin ip tan 9. Although 
these ratios can be expressed in terms of the spherical 
angles 9 and <p, we believe that presentation and visu- 
alization of the results using B' x and B' is simpler and 
more insightful than in the spherical angles [15l l25l |28| . 

The gauge phase © in Eq. © leads to interference be- 
tween interlayer tunneling amplitudes t c e l ^ r ' along the 
trajectory r(t). In Eq. y(t) oscillates with the period 
At = 2tt/oj c , whereas x(t) steadily increases, accumulat- 
ing the phase A<f> = edB y vpAt/Tic over one period. The 
average (e l ^'*')t vanishes unless A<j) = 2-rm, where n is 
an integer. This condition selects the Lebed magic an- 
gles B' y = n which in the spherical coordinates are 
sin<ys = n(b/d) cot 9 H3- Using Eqs. @, © and ©, we 
find the effective interlayer tunneling amplitude t c 



t c J n (B' x ) for 



where J„ is the Bessel function. 

AMRO result from a periodic modulation of the ef- 
fective interlayer coupling t c in Eq. @ due to interlayer 



quires that the flux of B y through the area, formed by the 
interlayer distance d and the electron trajectory period 
Air = vpAt, is n$0) where $0 = hc/e is the flux quan- 
tum. In addition, © oscillates as a function of BL 
with the period AB' X = it. These DKC oscillations jU 
are related to the flux of B x through the area bounded 
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B y =(B y / B z )(d/b) 



FIG. 2: The contour plot of angular oscillations in the nor- 
malized interlayer dc conductivity a c (B, 0)/o" c (0, 0), Eq. JHJ. 
The variables B' x and B' y are defined in Eq. (Sj. The lines 
with circles and squares indicate where interference between 
the two trajectories in Fig.Qa is constructive and destructive. 



by d and Ay = 4ti,c/ 'evpB zl the transverse width of 
the electron trajectory in Eq. More precisely, it is 
necessary to consider the distance between the turning 
points of an electron trajectory, as viewed along the vec- 
tor (B x ,B y ). This will be discussed in more detail from 
the momentum-space point of view. 

The interlayer ac conductivity <j c {uj) is given by a cor- 
relator of tunneling events at times t and t' [Til fi"o| 



<7 c (w) cx IZe t c 



j<f>(t')-i<fi(i) e (t'-i)(Uj-i/T) dt i 



(7) 

where r is a relaxation time. Substituting Eqs. and 
© in Eq. (JJ, we find 



tr c (0,0) 



E 

i— — 00 



Jn( B ' x ) 



1 + (uj c T) 2 (n - 



B' y T u/oJc) 2 ' 



(8) 



where a c (0,0) is the dc conductivity at B = 0, and 
the signs =F in the denominator originate from the ±vp 
sheets of the Fermi surface. Eq. © is in agreement with 
Refs. HE III El 

It can be applied to the mi- 
crowave measurements at lo ^ [35L l36j , but below we 
concentrate on the dc case u> = 0. When lu c t — ► 00, 
(6) only the term with n = B' y survives, and Eq. (JHJ reduces 
to cr c (B,0)/cr c (0,0) = (t c /t c ) 2 with t c from Eq. In 
Fig. we show the contour plot of a c (B, 0)/a c (0, 0) vs. 
B x and B' y calculated from Eq. JBJ for uj c t = \/50 w 7.1. 
The dc conductivity a c is maximal at the vertical stripes, 
labeled by the integer numbers n, which correspond to 
the Lebed magic angles B' y = n. Within the n-th vertical 
stripe, a c has alternating maxima and minima, indicated 
by circles and squares, which represent oscillations of J 2 
vs. B' x in Eqs. © and JHJ). Positions of these maxima 
and minima can be obtained from the Aharonov-Bohm 
interference in momentum space, as described below. 



B y =(B y / B z )(d/b) 



B x = ( B x /B z><W 



FIG. 3: Contour plot of the normalized interlayer dc conduc- 
tivity (J C (B, 0)/cr c (0, 0) vs. B' x and B' y calculated from Eq. 
1121 . Compared with Fig.|21 this plot takes into account addi- 
tional tunneling amplitudes t m along d + mb: t±i = t c /2 and 
t±2 — t c /4, which produce the Lebed oscillations at B x — 0. 



FIG. 4: Comparison between the normalized interlayer 
dc conductivities o c (B, 0)/u c (0, o) calculated from Eq. JHJ 
(dashed line) and Eq. 1121 (solid line) and plotted vs. B' x 
for By = 0. The DKC oscillations are reduced in the latter 
case, because of the additional tunneling amplitudes t m . 



Eqs. @ and ® show that, in the process of interlayer 
tunneling, the in-plane electron momentum changes by 



(qx,q v ) = (ed/c)(B y ,-B x ). 



(9) 



Thus, the Fermi surfaces of the two layers are displaced 
relative to each other by the vector q |l4Lll5| . as shown in 
Fig. Electrons can tunnel between the layers only at 
the intersection points fei, fc2, ^3, etc. of the two Fermi 
surfaces, where the conservation laws of energy and mo- 
mentum are satisfied. In the presence of B z , there is 
a phase difference between the two trajectories connect- 
ing the intersection points, which is proportional to the 
shaded momentum-space area Si > or S2 < in Fig. 
\T]p. The algebraic sum S± + S2 = q x (2Trh/b) depends only 
on q x = (ed/c)By. Constructive interference between fei 
and fc3 requires that (Si + S2)c/TieB z = 27rn, which is 
equivalent to the Lebed condition B' y = n. 

Interference between fci and &2 is controlled by the 
area Si. Introducing the dimensionless variable S[ — 
Sic/heB z , we find from Fig.^ 3 that 




^ ] w 

B' 



2 an -sin ( 



(10) 

Constructive interference requires that S[ = 27r(j + 1/4), 
where j is an integer, and the extra phase tt/2 appears 
because fci and &2 are the turning points on the Fermi 
surface, when viewed along the vector q. The lines with 
circles show where in Fig. El this condition is satisfied. 
Maxima of a c are achieved at the circled intersections of 
these lines and the integer vertical lines, where both Si 
and Si + S2 give constructive interference. These points 
correspond to the maxima of the Bessel functions in Eq. 
(JHJ. The lines with squares in Fig. show where the 
interference in Si is destructive (j is half-integer). At 
the intersections of these lines and the integer vertical 



lines, marked by squares, a c has minima, and the Bessel 
functions in Eq. (JSJ have zeros. There, a c — > at oj c t — > 
00, and resistivity p c = l/er c increases without saturation 
when B — > 00, whereas p c (B) saturates at the circles [32|. 
The maxima and minima of a c create a checkerboard 
pattern of oscillations 31] for \B' X \ > \B' y \ in Fig. HI 

The diagonal line BJ = B' in Fig.Elcorresponds to the 
third angular effect 0, HI HH • At this line, the points 
hi and &3 merge, and the area S2 shrinks to zero in Fig. 

For \B' X \ < \B' y \, the two Fermi surfaces do not cross 
in Fig. \\]p, so interlayer tunneling is suppressed, and a c 
does not show oscillations below the diagonal lines in Fig. 
121 However, this contradicts experiments [32],|38|], which 



show the Lebed oscillations of a r vs. B' at £L 



0. 



To improve the theory, let us consider a model with 
interlayer tunneling amplitudes t m between the chains 
shifted by m units in the y direction |3fj| . The tunneling 
displacement is d + mb, so the phase in Eq. J2Jl becomes 



h c 



(A z d + A y mb), 



Ay — B z x. 



(11) 



Comparing Eqs. Q and l|llfl. we see that results in this 
case can be obtained by substitution B y d — > B y d— B z mb 
and By — ► B' — to in the old results. Eq. © transforms 
into t m = t m J n ^ m (B' x ) for B' = n, and Eq. JSJ becomes 



<7c(0,0) 



E E i 

m n— — 00 



m ,J n — m 



+ M 2 ( Il -B'T^ c ) 2 



( 12 ) 

where = t 2 m j J2i tf- The contour plot of Eq. (|Ll)) can 
be obtained by shifting the plot in Fig. |3 by m units 
along the B' y axis and adding the shifted plots with the 
weights t'^. The resulting contour plot, calculated for 
to = t Cl t±i = t c /2, and t±2 = t c /A, is shown in Fig. 
IS At B x = 0, a c (B'y) has maxima for those directions 
B' = m where t m exists Oscillations of <r c vs. 
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B' x are smeared in Fig. [3J because the shifted maxima 
and minima of the checkerboard pattern in Fig. [21 add 
up out of phase. This is illustrated in Fig. 01 which 
shows that the DKC oscillations of <j c {B' x ) for B y = 
are much weaker for multiple t m . Moreover, a c does not 
have zeros at B' = n. Thus, when B — ► oo, p c (B) sat- 
urates on the integer lines in Fig. [3] but grows without 
saturation between the lines. Weak DKC oscillations of 
a c (B' x ) at By = and strong Lebed oscillations of cr c (B' y ) 
at Rr = correspond qualitatively to (TMTSF) 2 PF 6 
|23l 1371 138|. indicating that several t m are present. The 
opposite case, strong DKC and weak Lebed oscillations, 
is found in (TMTSF) 2 C10 4 0111111, suggesting that 
it has only one dominant t = t c |36| . The model pa- 
rameters can be determined by quantitative comparison 
between the calculated plots and experimental data for 
cr c (B' x , By). Fig. shows that the strength of the Lebed 
oscillations in a c vs. B' y increases when B ' ^ 0, in agree- 
ment with the Lee-Naughton experiment |27| . 

The amplitudes t m do not necessarily represent elec- 
tron overlap between distant chains. They may be effec- 
tive parameters in a model |39| . where e(k x ) has curva- 
ture, so vf depends on k x and varies along the quasiclas- 
sical trajectory Q. The resulting expression for a c has a 
form similar to Eq. Q12[l with some effective parameters 
t m , which themselves may depend on B pji^ . 

While Eq. may well describe the oscillatory part 
of er c , it often fails to describe the background, particu- 
larly in (TMTSF) 2 PF 6 HHHMHI, although there are 
variations with pressure and sample jijj. This remains 
one of the open problems, along with unusual tempera- 
ture dependence of resistivity 1371 and mysterious angular 
oscillations of the Nernst effect |41| . 

We presented a unified geometrical explanation of dif- 
ferent types of AMRO in Q1D conductors in terms of 
Aharonov-Bohm interference in interlayer electron tun- 
neling. We visualized a two-parameter pattern of oscilla- 
tions for generic magnetic field orientations using the nat- 
ural variables B' x and B' y . Quantitative comparison with 
experimental data plotted in this way is needed. This 
work was supported by the NSF Grant DMR-0137726. 
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